Magnetic Dispersion and Anisotropy in Multiferroic BiFeOa 
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We have determined the fuU magnetic dispersion relations of multiferroic BiFeOs. In particular, 
two excitation gaps originating from magnetic anisotropics have been clearly observed. The direct 
observation of the gaps enables us to accurately determine the Dzyaloshinskii-Moriya (DM) interac- 
tion and the single ion anisotropy. The DM interaction supports a sizable magneto-electric coupling 
in this compound. 

PACS numbers: 75.25.-j, 75.30.Ds, 75.50.Ee 



Multiferroic materials, in whieh spontaneous ferroelec- 
tric polarization and magnetic order coexist, have been 
investigated intensively not only due to their poten- 
tial industrial applications but also due to purely scien- 
tific interest about magneto-electric coupling in strongly 
correlated electron systems. For many geometrically 
frustrated magnets, ferroelectricity is mediated by the 
magneto-electric coupling. Several mechanisms have 
been proposed to explain the spin-driven ferroelectricity. 

a 

BiFeOs has a rhombohcdral structure (R3c) below 
~1100 K, where ferroelectricity appears. (This pa- 
per employs pseudo-cubic notation with a'^3.96 A and 
q;^89.4°.) The ferroelectricity is considered to primar- 
ily originate from displacements of the Bi'^"'" ions due to 
the lone 6s^ pair. Cycloidal magnetic order with the 
propagation vectors ti = ((5,-5,0), T2 = {S,0,—S), and 
'''3 = (0, — <^j4) develops below Tn ^ 640 K, as shown 
in Fig. [T] [J-Q The magnetic structure persists down 
to low temperatures, although the imcommcnsur ability 
S changes from -0.0045 at 5 K to -0.0037 at 600 K. 
[gHsI Because Tn is much higher than room temperature 
and because of the large spontaneous electronic polariza- 
tion (P — 100 /iC/cm^), @,[l3l this material has attracted 




many researchers and has been studied extensively. [11 1 
Although Tn is much lower than the ferroelectric Curie 
temperature in BiFeOs, several measurements show size- 
able magneto-electric coupling. For example, the mag- 
netic domain distribution can be controlled by applying 
an electric field. [J, [a, |l2| An abrupt decrease (up to 
~40 nC/cm^) in electric polarization was also observed 
in a magnetic field of about 20 T, [13l4l5| where a transi- 
tion from the incommensurate cycloidal structure to an 
almost commensurate structure with a weak ferromag- 
netic component is suggested, [la] These results suggest 
that the additional polarization below Tn is driven by 
the magneto-electric effect. However, it is important to 
clarify this mechanism from a microscopic point of view. 
Very recently, some inelastic neutron scattering studies 
have measured the spin- wave excitations using powder or 
single crystal sample. |17H19| Furthermore, the detailed 
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FIG. 1. (Color online) Magnetic structure of the Fe'^^ mo- 
ments in BiFeOs. The basic structure of coUinear G-type, in 
which nearest-neighbor spins align antiferromagnetically, is 
shown. The dotted lines correspond to the 2x2x1 unit cells 
of the rhombohcdral structure (_R3c). The pseudo cubic unit 
cell that is used in this paper is shown by the thin solid lines. 
The spin structure is spiral along the [1,-1,0] direction with 
a long period of ~640 A. The nearest-neighbor (Ji) and next- 
nearest-neighbor interactions (J2) as well as the directions of 
the DM vector (D) and the spontaneous electric polarization 
(P) are also shown. 



spin Hamiltonian including the magnetic anisotropy was 
discussed to explain magneto-electric coupling in this 
compound by Sosnowska et al. [20| and Jeong et al. [l9| . 
However, a detailed analysis based on the direct obser- 
vation of the magnetic anisotropy, which is necessary to 
discuss the magneto-electric coupling, has not yet been 
performed. 

We performed inelastic neutron scattering experiments 
on a single crystal of BiFeOs. We have determined the 
full magnetic dispersion relations of the spin- wave excita- 
tions in this compound. In particular, low-energy gapped 
excitations have been detected for the first time by high 
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FIG. 2. (Color online) Typical constant-Q and constant- 
energy scans along [h, h, h], [h, h, h], [h, h, 0], and [0, 0, h] mea- 
sured on HB-1 using thermal neutrons at T—200 and 300 K 
in BiFeOa. The solid curves are the results of fits of a convo- 
lution of the resolution function with Lorentzians. 



energy resolution experiments. The direct observation of 
the excitation gaps makes it possible to determine accu- 
rately the magnetic anisotropics due to Dzyaloshinskii- 
Moriya (DM) interaction and single ion anisotropy. Our 
detailed analysis has revealed that the coupling constants 
are Ji=6.48 meV, J2=0.29 meV, i:)=0.1623 meV, and 
i<ir=0.0068 meV, where Ji, Ji, D, and K are nearest- 
neighbor, next-nearest-neighbor, DM interactions, and 
single ion anisotropy, respectively, which reproduce the 
cycloidal spin structure. We also measured the temper- 
ature dependence of the dispersion relations below room 
temperature. Although no drastic change was observed, 
the spin- wave are slightly softened above 200 K. 

A single crystal of BiFeOs was grown using the travel- 
ing solvent floating zone (TSFZ) method by laser heat- 
ing, as described in Ref. |2l|. The dimensions of the sin- 
gle crystal is ~5(/)x40 mm^. The effective mosaic of the 
single crystal is about 0.8° with the spectrometer config- 
urations described below. The inelastic neutron scatter- 
ing experiments were carried out on the thermal triple- 
axis neutron spectrometer HB-1 and the cold triple- 
axis neutron spectrometer CTAX, installed at HFIR at 
ORNL. Neutrons with a final energy of 14.7 meV and 
3.5 meV were used, together with a horizontal collimator 
sequence of 48'-80'-S-80'-120' and guide-open-S-80'- 
open on HB-1 and CTAX, respectively. Contamination 
from higher-order beams was effectively eliminated us- 
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FIG. 3. (Color online) Low-energy magnetic excitations at 
the magnetic zone center (1/2, —1/2, 1/2) measured at T=5, 
100, 200, and 300 K on CTAX using cold neutrons. The 
solid curves are guides to the eye. The inset is a schematic 
figure to show the relation between the dispersions and the 
instrumental resolution (shaded ellipsoid). 



ing PC and Be filters on HB-1 and CTAX, respectively. 
In both experiments, the single crystal was oriented in 
the {HHO)-{OOH) scattering plane and was mounted in 
a closed-cycle ^He gas refrigerator. 

Figure [2] shows the typical inelastic neutron spectra 
along [h,h,h], [h,h,h], [h,h,0], and [0,0, /i] in BiFeOa 
measured on HB-1. Since the dispersion is steep and 
the instrument resolution in Q is not sufficient to re- 
solve the dispersions arising from the incommensurate 
positions, the observed spectra are considered to be al- 
most the same as those expected for the commensurate 
G-type structure, as shown in the inset of Fig. |3l and it 
is difficult to resolve the two peaks at -\-q and ~q below 
_E'~20 meV in constant-energy scans. In order to com- 
plete the spin-wave dispersion relations, we mostly used 
constant-energy scans between ^^20 and '-^60 meV and 
constant-Q scans between ^10 and ~20 meV and around 
the magnetic zone boundary. Below 10 meV we used the 
constant-energy scans measured on CTAX. In order to 
determine the peak positions of the spin- wave excitations 
in the energy-Q space, the CTAX data were fitted using 
the Gaussian function without convoluting with the res- 
olution function. The HB-1 data were fitted using the 
Lorentzian function A/[{E — _Eo)^+r^] with r=:1.5 mcV, 
where A and Eq are constant and peak position in energy, 
respectively, convoluted with the instrumental resolution 
function. As shown in Fig. [51 the model function repro- 
duces the observed spectra reasonably well. Therefore, 
the broad peak widths in the constant-Q scans are not 
intrinsic but primarily originate from the steep dispersion 
around the zone center (<^30 meV) and from the insuf- 
ficient instrumental energy resolution at higher transfer 
energies (>~30 meV). 

The experiments were carried out primarily at T=:200 



80 



[h,-h,0] 
T=200 K 




[h,-h,h] 
T=200 K 

- (c) 

T 

• / 

- 7 

/ 




- 


/ 


• [h,-h,h] 
[h,-h,-h] 




, , , , 1 , , , , 1 , , , 





1 



0.1 0.2 0.3 0.4 0.5 0.05 0.1 0.15 0.2 0.25 

h (r.l.u.) h (r.l.u.) 




0.1 0.2 0.3 0.4 0.5 
h (r.l.u.) 



0.05 0.1 0.15 0.2 0.25 
h (r.l.u.) 



FIG. 4. (Color online) Magnetic dispersion relations along 
[/i, /i,0] (a), [0,0, ft] (b), [ft, ft, ft] (c), and [ft, ft, 2ft] (d) in 
BiFeOs. The solid curves are spin- wave dispersion relations 
calculated with Ji=6.48 meV and J2=0.29 meV. 



K. Temperature dependence of the dispersion along 
[0,0, ft.] was also measured below 300 K. As shown in 
Figs. [2le),[2][f) andlU^b), the dispersion becomes slightly 
softened at 300 K with a change of ^2 meV at the zone 
boundary. The change of dispersion is almost negligible 
below 200 K (not shown). 

Figure [3] shows the temperature dependence of the low 
energy magnetic excitations at the magnetic zone cen- 
ter (1/2, -1/2, 1/2) measured on CTAX. Because there 
are several magnetic domains which have different prop- 
agation vectors, there are 12 incommensurate magnetic 
Bragg positions around (1/2, -1/2, 1/2), i.e. (1/2±(S, 
-1/2t<5, 1/2), (l/2±(5, -1/2, 1/2t^), (1/2, -l/2±5, 
1/2t5), (1/2±<5, -1/2±5, 1/2), (1/2±J, -1/2, 1/2±(S), 
and (1/2, -l/2±(5, l/2±(5). Although only the excita- 
tions arising from (l/2±5, —1/2^5, 1/2) arc observed at 
higher energies, some of the excitations are superposed 
around the zone center due to the instrument resolution 
in Q. Even so, gap energies can be determined because 
they are the same for all domains. At 300 K, we found 
a sharp peak at ~1.1 meV and a broad peak at ^2.5 
meV. Those excitation peaks are considered to be mag- 
netic in origin, because they are dispersive and contin- 
uously connected to the spin-wave excitations described 
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FIG. 5. (Color online) The relation between D and K that 
reproduce the cycloidal magnetic structure with a long period 
of ~640 A. The gap energy at the zone center due to the DM 
interaction Adm is also plotted as a function of K. 



above. Furthermore, (1/2, —1/2, 1/2) corresponds to a 
zone boundary of the chemical reciprocal lattice unit so 
that low-energy phonons are not expected to be observed 
at this position. Those peaks are broader than the in- 
strumental energy resolution ~0.2 meV around i5=l~2 
meV. Since the dispersion curve is steep, a slight tail is 
observed above the gap energy, as shown in the inset of 
Fig. [3l The gap energy of the lower excitation mode is 
estimated to be ~ 1.1 meV, whereas that of the higher ex- 
citation mode is estimated to be ~2.5 meV, which is dif- 
ficult to determine more accurately. As will be described 
below, these gaps probably originate from the DM in- 
teraction and the single-ion anisotropy that give rise to 
an easy-plane anisotropy in the plane defined by [1, 1, 1] 
and the spiral direction [1,-1,0] with a finite in-plane 
anisotropy along [1,1,1]. The scattering intensity below 
0.8 meV is considered to originate from the magnetic ex- 
citations from the almost gapless mode. With decreasing 
temperature, the scattering intensity decreases, follow- 
ing the Bose factor. The lower gap energy just slightly 
increases, which is consistent with the result that the 
incommensurability does not change considerably below 
300 K. [6|-|8| It is noted that the magnetic anisotropy is 
small, as expected from the absence of the orbital degree 
of freedom in the Ye'^ ions (3d^ 5=5/2). 

Figure 2] shows the full magnetic dispersion relations of 
the spin- wave excitations along [/i, h, 0], [0, 0, /i], [/i, /i, /i], 
and [/i, /i,2/i]. The maximum energy of the spin- wave 
excitation is ~70 meV. In order to analyze the spin-wave 
dispersion relations observed, we assumed the following 
effective spin Hamiltonian. [19|, |20|, |22| 



ij = Ji ^ s., • s, + ,h J2 ^^ ■ ^1 

n.n. 7i.ri.7i. 

[1.-1,0] i 



The first and the second terms represent exchange in- 
teractions between nearest-neighbor and next-nearest- 
neighbor spins, respectively. The third and the fourth 
terms originate from the DM interaction and the single- 
ion anisotropy, respectively. A moment size of 4.1^^ at 
200 K was used in the calculations, since the saturated 
moment was reported to be 4.34/iB and it is reduced by 
^5% at 200 K.[23[ Calculations for the spin dynamics 
were performed by including both the DM interaction D 
and anisotropy K in the Heisenberg Hamitlonian. For a 
fixed period of the spin helix, K depends on D. Since 
5 ~ 0.0045 is approximately equal to 1/222, a unit cell 
of length 222 was used to evaluate the anharmonic con- 
tributions to the spin helix, which was expanded in odd 
harmonics of a fundamental wavevector Q. [2J| The spin 
excitations were evaluated by performing a 1/ S expan- 
sion in the rotated frame of reference for each spin in 
the unit cell. The equations-of-motion for the spin op- 
erators were solved by diagonalizing a 444-dimensional 
matrix. This technique is described in more detail in 
Ref. [25|. Zone folding then generated the two gap fre- 
quencies, which depend on the anisotropy K. 

We first consider the overall dispersion relations, which 
are mostly caused by the first and the second terms. 
The solid curves are spin- wave dispersion relations calcu- 
lated with Ji=6.48 meV and J2=0.29 meV. These values, 
which describe the observed dispersions reasonably well, 
are consistent with Ji=4.38 meV and J2=0.15 meV from 
Ref. [l9|, where a moment of 2hb\/S{S + 1) = 5.8fiB 
or 1.44 times the T = 200 K moment was used. Our 
results predict the band maximum at 72.8 meV, which 
is also consistent with 72.5 meV in Ref. [19|. Magnon 
density of states measured with polycrystalline sample 
is supposed to show a peak around the band maximum. 
[26| The peak position was experimentally determined at 
68.2 meV at 20 K [l3| and at 65 meV at 300 K [ll|. These 
are slightly lower than the predicted values probably be- 
cause the instrumental resolution broadens and lowers 
the asymmetric peak from the van Hove singularity. 

Magnetic anisotropics should be included to explain 
the excitation gaps around the magnetic zone center, al- 
though they do not affect the excitations above ~5 meV. 
The DM vector is assumed to point along the [1,1,-2] di- 
rection and the summation was done for spin pairs along 
the [1,-1,0] direction. It was reported that there is a fi- 
nite easy-axis anisotropy along the [1, 1, 1] direction due 
to the single-ion anisotropy. [a, [24I These terms give 
rise to excitation gaps around the magnetic zone cen- 
ter shown in Fig. [31 D and K cannot be chosen in- 
dependently to reproduce the cycloidal magnetic struc- 
ture with a long period of ~640 A previously observed. 
[3|, |lO[ Furthermore, the single-ion anisotropy lifts the 
degeneracy of the excited state. The relation between D 
and K is plotted in Fig. [51 In the region of A'>0.0222, 
the collinear G-type magnetic structure becomes stable. 
Since the lower excitation gap, which is sharp and well- 



defined, is observed at 1.10±0.05 meV, K is estimated 
to be 0.0068±0.0007 meV. Then, D is estimated to be 
0.1623±0.0022 meV. Consequently, the higher excitation 
is predicted to be at 2.33 meV, which is comparable to 
the experimental result, although the exact peak posi- 
tion is difficult to locate experimentally. The gap energy 
due to the easy-axis anisotropy K is predicted to be less 
than 0.1 meV. Therefore, it is difficult to observe the 
anisotropy gap experimentally. The magnetic excitation 
observed below 0.8 meV corresponds to this spin- wave 
mode. The lower gap energy is considered to change lin- 
early with the Fe moment, which decreases by ~5% from 
100 to 300 K. mill This is consistent with the behavior 
that the lower gap energy slightly decreases by ^^3% from 
100 to 300 K, as shown in Fig. [3 

The K term gives rise to an anharmonicity of the 
cycloidal structure, which was previously suggested by 
NMR and neutron diffraction studies, [61, [sj, |27[ Higher 
harmonics of the incommensurate magnetic peaks are in- 
duced due to the anharmonicity. From /ir=0.0068 meV, 
the ratio of the intensities of the first to the third harmon- 
ics is estimated to be /i//3~120, which is between ~500 
at T=5 K l6| observed in the neutron diffraction mea- 
surement and ~25 at T=5 K estimated from the NMR 
measurement. [27[ Further studies are required to under- 
stand the discrepancy between these measurements. 

It is interesting that the zone-center gaps around 1 and 
2.5 meV are consistent with the magnon modes observed 
in Raman scattering and terahertz spectroscopy measure- 
ments, [22l . I28l43l| in which the modes are considered to 
be electromagnons. This indicates that the zone-center 
spin- wave excitations were observed previously, although 
the electromagnon modes are observable at q=0 in prin- 
ciple. Our results give important information to under- 
stand the coupling in those measurements. 

The direct observation of the excitation gaps enables 
us to determine an accurate value for the DM interaction. 
The presence of the DM interaction suggests that the ad- 
ditional ferroelectric polarization below Tn is caused by 
the magneto-electric effect predicted by the spin current 
mechanism. [S^l Therefore, the cycloidal spin structure 
in BiFeOa is not caused by geometrical magnetic frustra- 
tion as in many other multiferroic materials but rather 
by the relatively large DM interaction. Hence, pertur- 
bations that increases the DM term should enhance the 
spontaneous electric polarization. 

In summary, from the direct observation of the excita- 
tion gaps, we determined all the magnetic coupling pa- 
rameters, including the DM interaction and the single-ion 
anisotropy, which are consistent with the observed anhar- 
monicity of the cycloidal magnetic structure in multifer- 
roic BiFeOa. Based on the DM interaction, we conclude 
that the additional ferroelectric polarization below Tn is 
caused by the sizable magneto-electric effect. 
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